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Clusters of carbon atoms have been studied with large-scale ab initio calculations. Planar, 
single-sheet graphite fragments with 6 - 54 atoms were investigated, as well as the spherical C 6 o 
"Buckminsterfullerene’ molecule. Polycyclic aromatic hydrocarbons (PAHs) have also been 
considered. Thermodynamic differences between diamond- and graphite-like grains have been 
studied in particular. Saturation of the peripheral bonds with hydrogen is found to provide a smooth 
and uniform convergence of the properties with increasing cluster size. For the graphite-like 
clusters the convergence to bulk values is much slower than for the three-dimensional complexes. 


Introduction. 


Due to their scientific and technical importance, carbon clusters have long been subject to a 
variety of experimental[1-24] and theoretical[25-5lj investigations. The areas of interest are of 
great diversity, including heterogeneous catalysis, the composition of graphite vapor, soot for- 
mation during combustion[15] and early stages of synthetic natural gas production from coal. 
Several small clusters have been identified spectroscopically in the atmospheres of red giant stars 

and in comet tails[6-9]. ... 

The experimental work reported for carbon clusters ranges from the early investigations of 
Honig and Drowart et al on small clusters[1 -3] to the recent studies in the groups of Smalley and 
Kaldor[10-16] on larger clusters with up to 200 atoms. 

A striking observation that lacks a satisfactory explanation is the existence of 'magic numbers', 
i.e. the fact that in a distribution of clusters some species with a certain number of carbon atoms 
are much more abundant than others. The exact clustering mechanisms are not completely 
understood, and, as noted by e.g. Rohlfing et al.[10], the origin of the observed distribution of clus- 
ters may depend upon instrumental factors. Accounting for this fact, however, there still seems to 
be a preference for clusters with certain numbers of atoms which cannot be explained solely as due 
to the experimental conditions. 

Evidently, several aspects of this exciting area are difficult to study with experimental tech- 
niques. The different species are short-lived, reactive, and exist only under rather extreme condi- 
tions. These are conditions under which theoretical studies can contribute a lot to our under- 
standing. Theoretical work has indeed been reported on smaller clusters with n=2-10[42-50] as 
well as on some of the larger ones[25-41]. The present work reports ab initio calculations for a 
number of large carbon clusters of relevance for the chemical problems addressed above. 
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Ab initio investigations on systems of the size considered here are usually not considered a 
routine procedure. Calculations which account for electron correlation in a meaningful fashion are 
still beyond reach for most of these molecules. With recent progress in hardware and software de- 
velopment it is now technically and economically feasible to study quite large molecules using the 
Hartree-Fock method. That approach has also the advantage of being size consistent, which is a ne- 
cessity for this type of study when results for molecules of vastly different size are to be compared. 
This allows for the selection of study species which may at first be relatively small, but which may 
be gradually increased to sizes that eventually strain the capacity of the computing equipment. 
Size-consistency allows for an extrapolation of these results to even larger molecules, eventually 
to bulk systems, which are far beyond reach for a direct application of these computational methods. 


Graphite- like Fragments 


Due to the enormous number of possible structures for all but the few smallest clusters, quite 
restrictive assumptions about their geometry have to be made. For those clusters where well-de- 
fined equilibrium structures do exist, these are likely to possess a non-trivial point group sym- 
metry (in many cases the highest possible symmetry). It there fore seemed justified to focus the 
study on high-symmetric systems. Symmetry can also be used to simplify the calculation of elec- 
tronic structure, and reduces the number of geometrical degrees of freedom to be optimized. In or- 
der to further confine the study, only graphite-like systems consisting of pericondensed aromatic 
six-rings are considered. 



Figure 1 . Graphite-like fragments of D 6h symmetry. C 6 (1), C 24 (2) and C 54 (3). 


In an initial series of calculations, the systems C 6 , C 24 and C 54 (D 6h symmetry, 1, 2 and 3 in 
Fig. 1) were computed in a high spin form with all unpaired electrons having parallel spins (Total 
spin S=3,6, and 9). All carbon-carbon bond distances were kept equal at 1.42 A, the distance ex- 
perimentally found in bulk graphite(55). The total energies computed, -226.5321 Hartree for C 6 , 
-907.0085 Hartree for C 24 and -2041.5720 Hartree for C 54 , can be interpreted in terms of co- 
hesive energies of 71.0, 97.1 and 106.5 kcal/mol per carbon atom, respectively. 
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ADDarentlv the cohesive energy of these clusters shows a very slow convergence with the size 
of the P molecule. This should not be surprising, since the number of unsaturated valences dangling 

honds" Der carbon atom is one in 1, one-half in 2 and one-third in 3. , , . 

Assuming an approximately constant energy per C-C bond, the above trend is understandable. 

With clusters on the above general type, the number of carbon atoms is 6N , the d *"' 

alino bonds is 6N and the number of C-C bonds is 9N 2 -3N. The energy per bond shows a smoother 
S.^hf numbers big 71.0. 77.6 and 79.9 kcal/mol. respectively. Alternatively, the energies 

can be fitted to a two-parameter expression of the form 


E tot = 6N 2 E C + 6NE db 


(D 


where Ec(energy per carbon atom) and E db (energy per dangling bond) are adjustable parameters. 
The optimum coefficients are E c =-37.8366 Hartree and E db = +0.0891 Hartree. For large values 
of N the total energy per carbon atom will converge towards E c , and accordingly the cohesive energy 
per carton atom in graphite can be expressed as the difference between E c and the atom, energy of 
Srhnn r ^7 6*171 Hartree) leading to a value of 0.1995 Hartree or 125 kcal/mol. (i.e. 83.3 
kcSmol ce- C C tend since in g?aphite there are 3/2 bonds per carbon atom). The Parameter 
E for Which a value of +0.0891 Hartree is derived, can be interpreted as the energy loss (note 
^positive sign of E db ) for introducing a dangling bond in a graphite-like sheet of carbon atoms. 
Since breaking a C-C bond would result in two dangling bonds, the value of E db corresponds to a loss 
of bond energy of 112 kcal/mol per bond, somewhat higher than the bond energy derived from E c . 
One miaht have expected the electrons unpaired by the bond-breaking to redistribute into o er 
bonding 9 orbitals orjn valence-bond language, that the broken peripheral bonds would rehybridize 
to minl°ze tffe destabilization of the framework. The opposite turns out to be true, though 
breaking some bonds weakens others due to an interrupted electron delocalization, a phenomenon 
SKSSd by the fact that the C-C bond in cyclic C 6 is considerably weaker than m graphite. 

For the larger clusters C 24 and C 54 all carbon atoms are not symmetry equivalent and fluc- 
tuations in the local properties from atom to atom may occur. It is of some interest to i study 'the en- 
tails of that fluctuation. For instance, it is important to know how large the molecule needs to be 
before its interior is reasonably close to that of a single-sheet graphite. The population analysis for 
C 24 shown in Fig. 2 demonstrates that the charge originating from the unpaired electrons is located 

\o \ large extent in the peripheral carbon sigma-bonding orbitals. 


0.969 



Figure 2. Spin-density populations for C24 in a high-spin form (S-6). 

So far only high-spin systems have been considered. However, when neighboring dangling 
bonds interact, spin-pairing may occur and lead to the formation of new bonds, and hence to a lower 


247 


energy. In the case of C 24 all adjacent dangling bond electrons can be paired up, which would result 
in a singlet state of the molecule. A partial structure optimization on this state leads to very short 
peripheral carbon-carbon bonds (1.20 A) with essentially triple bond character. All other bond 
lengths change very little despite the strain introduced by shortening one bond of the outer six- 
membered ring by more than 0.2 A. For C 54 a similar structure may be obtained. Pairing up elec- 
trons from adjacent carbon atoms of the same six-membered ring, and leaving all the other 
(isolated) dangling bond electrons with parallel spins, a low spin state (S=3) would be obtained. 
This mechanism is supported by an actual calculation on the low spin form, in which the spin den- 
sity is almost entirely localized on the carbon atom that would carry the dangling bond in a simple 
valence orbital picture, (see Fig 3). 

low-spin (S=5 ) 



Figure 3. Spin density distribution for C 54 in a low-spin form (S=3). 

Energy calculations show that the low spin state of is indeed about 1 .3 eV lower in energy 
than the high spin state, at least for the molecular geometries used here, whereas for C 24 nearly 
identical energies for the two states are found. Inclusion of electron correlation would favor the low 
spin form further, though hardly enough to significantly affect a comparison of cohesive energies 
between various forms of carbon. 

The ionization potentials (Koopmans' values) for the high spin and the low spin form of C 54 are 
6.8 and 7.3 eV respectively, values that would both fit into the bracket of IP’s reported for carbon 
clusters with 40 to 100 atomsjll]. The corresponding values for C 24 are higher, 8.0 and 8.3 eV 
respectively. 


Polvcvcll c Aromatic Hydrocarbons. 


The above studies show that properties in the interior of the graphite-like clusters do not 
converge readily to bulk values. Even for the C 54 cluster, the populations in the center show a sig- 
nificant deviation from neutrality. Despite the possibility of exploiting the high point-group sym- 
metry of these systems, the complicated open-shell structure places severe limits on the maximum 
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size of the clusters that can be studied. Therefore, the calculation of a graphite cluster large enough 
to ensure that the calculated properties have converged to bulk values would be computationally 
unfeasible. 

If the focus of interest is on the carbon clusters themselves, then of course no substitute sys- 
tem can be used. However, for studying the convergence of properties one can minimize the termi- 
nation effects by saturating the dangling bonds in the simplest possible way, i.e. with hydrogen. By 
that approach one can both avoid the problem of handling an excessive number of open shells, and 
obtain a series of molecules that converge towards bulk properties more smoothly than the bare 
carbon clusters. 



Figure 4. The PAHs C 6 H 6 (4), C 24 Hi2(5), C 54 H 18 (6). C 96 H 24 (7), and C 15 oH 30 (8) 

The series of PAHs shown in Fig. 4 have been studied in the present work, and should allow for 
an extrapolation towards an infinite single layer of graphite. The total energies of these systems 
with the general formula C 6 n2H 6 n are presented in Table 1. For all molecules except 8 an overall 
C-C bond distance was optimized. This seems to converge smoothly towards 1 .406 A, a value some- 
what shorter than that experimental distance of 1.421 A found in bulk graphite[52]. The Hartree- 
Fock approximation is the main source of this discrepancy. It should be noted, however, that the 
bond distance in a graphite mono-layer is expected to be somewhat shorter than the bulk value[53]. 

The total energy of these molecules can be accurately fitted with a two-parameter expression 
similar to that given in Equation 1 : 


E tot,est - GN^Eq + 6NEh 


(2) 


Table 1 shows that the total energies of the systems studied here are reproduced by Eq. 2 with 

an error smaller than 0.2 kcal/mol per carbon atom. 

The optimum coefficients in Equation 2 are Ec=-37.8363 and Eh=- 0.5685 Hartree. Eq is 

very close to the value -37.8366 found for the pure graphite clusters. 
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Figure 5. HOMO and LUMO orbital energies for C 6N 2H 6N , plotted vs 1/N. 

The value of E H corresponds to a contribution of 47.1 kcal/mol to the total energy of C 6N 2H 6N for 
each bonded H-atom (the hydrogen atom has an energy of -0.4935 Hartree in the basis set used). 



Alternatively, the parameter values can be interpreted as 83.3 kcal/mol per C-C bond, 88.7 
kcal/mol per C-H bond. 

The value 1 25 kcal/mol represents an upper bound to the cohesive energy per carbon atom in 
graphite, since the interaction between layers in the bulk has not been accounted for. Given the 
relatively large distance and the physical properties of graphite, the inter-layer interaction en- 
ergy is estimated to be < 5 kcal/mol. 

This I 03 V 0 S us with 3 computod rosult in loss thsn sstisfsctory sgroomont with tho oxpon- 
mental value of about 170 kcal/mol[54]. The neglect of electron correlation and the limited basis 
set used are the most important sources of the discrepancy. In a previous study on mono-layer 
graphite[53], basis set effects were found to lead to a significant underestimation of the cohesive 

onorgy. 

The orbital energies show a surprisingly slow convergence as the size of the molecule in- 
creases. Even for the largest system studied, the HOMO-LUMO gap is about 4 eV. From the values 
shown in Fig. 5 one may expect semiconducting behavior to occur at N=8-10, i.e. for clusters with 
400-600 atoms. At that point a conventional Hartree-Fock approach would of course break down. 

The HOMO and LUMO orbital energies seem to converge towards a value of about -4 eV, which 
would be an estimate of the work function for single-sheet graphite. This is in reasonable agree- 
ment with the work function of 4.9 eV experimentally found for bulk graphite[52]. 


Dlamond-l ike Clusters. 


From the two-dimensional, graphite-like clusters, the extension to three-dimensional struc- 
tures is obvious. Symmetric structures developed in a similar fashion to the planar systems would 
grow in three dimensions with increasing N, and the number of atoms would increase faster. In this 
work clusters of T,j symmetry were studied, resembling a small fragment of a diamond structure. 
Only systems with saturated external bonds were considered. The number of carbon and hydrogen 
atoms in such a structure is given by 

Nq = N(4N 2 -1)/3 
N h = 4N 2 


Here we present some results from studies on systems with N-2,3, and 4. (N-1 would cor- 
respond to the methane molecule). The molecules are shown in Fig. 6 

For all three molecules, an overall C-C distance was optimized. The values obtained for the 
three systems are 1.542A, 1.540A and 1.539A, smoothly converging to a value slightly shorter 
than the equilibrium distance in diamond. 

Fitting the energies to an expression similar to Eq 2 gives E c - - 37.8360, Eh - -0.5732 for 
the tertiary hydrogens, and E H - -0.5742 for methylene carbons. The extrapolated cohesive energy 
is 125 kcal/mol, the same as for graphite. A breakdown of the bond energies per bond gives a C-C 
bond energy of 62 kcal/mol. The C-H bond energies are 82 and 81 kcal/mol for the methylene and 

the tertiary hydrogens, respectively. _ . . .. . . . . ,. 

The convergence pattern for the orbital energies is quite different from that found in the 
graphite-like clusters as shown in Fig. 7. The extrapolated band gap is somewhat uncertain but 
seems to lie around 10 eV, and is definitely larger than the experimental value of about 4.6 eV. This 
difference is mainly due to neglect of electron correlation in our calculations. 
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Figure 7. HOMO and LUMO energies for the diamond-like clusters, plotted vs 1/N. 


Another simple form of elemental carbon would be chains formed from carbon atoms. As a 
prototype model a single-stranded chain is most suitable. If branching were to be considered L all 
intermediate forms up to and including the diamond- and graphite-like clusters wou d be mckided. 
For non-branched chains, the two variants to choose from are a system of alternating singly and 
triply bonded carbon atoms (poly-ynes), and a system with all double bonds (cumulenes). Cumu- 
lene V structures are assumed to be the preferred ones for odd-membered chains, whereas the even 
ones may have some poly-yne character. Recent studies on linear Cq show that a cumutene-like 

structure is preferred, both at the SCF level and when correlation is accounted for[50J. 

As an alternative to linear chains, a simple ring-shaped molecule should also be considered. 
Reaardless of whether poly-ynes or cumulenes are considered, the gain in bond energy due to the 
new bond formed should be relatively independent of the size of the chain. In contrast, the strain 
involved in forming a ring is inversely proportional to the number of atoms, assuming a harmonic 
C C-C bending potential. Clearly, ring structures should be thermodynamically preferred for 
chains above a certain size (even though their formation might be kmetically or statistically un- 
favorable). In fact, already for as small a system as C 6 the ground state structure is found to be 
cyclic, although the omission of electron correlation would predict a linear, cumulene-like geom- 
etrv to be marginally lower in energy[50]. 

We have chosen to study the C 20 poly-yne ring assuming it to be a reasonable representative 
for a large chain molecule. The bond lengths optimized at the Hartree-Fock level are 1 .37 A and 
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1 .20 A, respectively. The same values are found for the experimental bond lengths in 1 ,3-buta-di- 
yne CH 3 -C=C-C=C-H, indicating that there is very little n conjugation in the system. 

The cohesive energy per carbon atom in a poly-yne ring is only 99.1 kcal/mol, clearly lower 
than the value in e.g. planar C54. In addition to the unfavorable energy one would anticipate a long 
and complicated route of formation for a chain-like molecule when starting from graphite, and it 
does not seem likely that any of the larger clusters observed experimentally would have a linear or 
cyclic chain structure. 


Non-planar Clusters Containing Flve-membered Rinas. 


The existence of a stable Cgo molecule, Buckminsterfullerene was recently proposed[14). In 

the suggested structure, twelve regular pentagons and twenty hexagons are connected to form an al- 
most perfectly spherical pattern of icosahedral symmetry. 

Previous theoretical work seems to support the existence of this amazing molecule. Several 
calculations[25-41] on Buckminsterfullerene at various degree of approximation find a reso- 
nance stabilization larger than in benzene, and only slightly smaller than in graphite. From these 
calculations it has been concluded that the molecule is likely to exist. A recent MNDO study [41] 
concludes that large polyhedral complexes (> 40 atoms) are more stable than the corresponding 
graphite-like ones. However, the predictive power of these approximate methods is limited, in 
particular for problems where planar and non-planar systems are compared. The interpretation of 
the original experiment postulating the existence of Buckminsterfullerene has recently been chal- 
lenged[1 1 ], and ab initio calculations of the relative energies involved could be used to resolve some 
of the controversy. 

The optimized geometry has twelve regular five-membered rings with an edge of 1.453A. Each 
of the thirty hexagons shares three sides with pentagons, whereas the three others, shared with 
adjacent hexagons, are substantially shorter (1 .369A). Our calculated HOMO-LUMO splittings are 
around 8 eV, both at the minimal-basis and double-zeta level. The results of our calculations on 
Buckminsterfullerene were reported fully in Ref 39, and are summarized in Table 3. 

The ionization potential (7.9 eV) falls right outside the bracket of experimental IP's reported 
for carbon clusters with 40 to 100 atoms (6.42 eV SIP £ 7.87 eV[11]). Inclusion of correlation 
effects will lower the calculated ASCF IP by 0.25 to 0.50 eV, so that the corrected IP will be at the 
upper end of the experimental IP-bracket. Due to the diffuseness of the n orbital from which an 
electron is removed, the correlation error in the ASCF value will be smaller than in cases where an 
electron is removed from a well localized bond. In those cases a correction of 1 eV is usually 
applied. 

The cohesive energy per carbon atom in Buckminsterfullerene, 114 kcal/mol, is 11 kcal/mol 
lower than the value extrapolated for graphite, but 7 kcal/mol higher than the value computed for 
the C 54 graphite-like sheet or 14.5 kcal/mol higher than for the C20 poly-acetylene ring. Accord- 
ingly, the heat of formation for Buckminsterfullerene is 650 kcal/mol, suggesting a high degree of 
instability with respect to graphite. 

At the Hartree-Fock level, comparisons between systems with very different molecular and 
electronic structures have to be made with some care. Strained geometries of the type found in C 60 
are usually difficult to describe without polarization functions. The basis set used is therefore 
expected to slightly favor the planar systems. Repeating the calculations with d-functions added 
would be unnecessarily time-consuming, and another approach has instead been used. The bowl- 
shaped corannulene molecule (C^H^) has a geometry resembling a fragment of a C 60 ball, and the 
strain in the C-framework ought to be similar as well. 

To investigate the basis set effect in connection with geometric strain in Buckminsterfullerene, 
we have performed calculations with and without d-type functions on corannulene. The d-functions 
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improve the binding with about 10 kcal/mol per C-atom in C^H^; the corresponding improvement 
in coronene (C 24 H 12 ) is 8 kcal/mol. Within reasonable limits, these numbers are not crucially 
dependent on the d-exponents used(0.4-0.6). One can therefore estimate that the inclusion of d- 
orbitals in Buckminsterfullerene would decrease the calculated heat of formation by about 120 
kcal/mol. 

Even though the predicted relative instability of Buckminsterfullerene would be somewhat di- 
minished by such an extension of the basis, the main conclusion remains; the molecule is still 
about 530 kcal/mol higher in energy than a corresponding fragment of bulk graphite. 

Another conclusion that can be drawn from the calculations on corannulene is an estimate of the 
loss of energy due to incorporation of a five-membered ring into a structure of condensed six- 
rings. Eq. (2) predicts a total energy of -762.4110 Hartree. for C 20 H 10 , the energy actually com- 
puted is -762.2656 Hartree., and the loss due to a five-membered ring is thus about 0.14 
Hartree. This is significantly larger than loss of 0.09 Hartree. predicted from Eq. (1) due to the 
occurrence of a dangling bond in the structure. It can therefore be concluded that the reducing the 
number of dangling bonds in 2 and 3 above by introducing five-membered rings on the periphery 
is not likely to lead to more stable structures. 

In conclusion, these calculations suggest that the Buckminsterfullerene molecule is thermo- 
dynamically unstable with respect to graphite by about 23 eV. In general, thermodynamically un- 
stable molecules can often be readily synthesized, and frequently exist with virtually infinite 
lifetimes. The unusually large instability observed in this case might seem almost prohibitive, 
however no other good candidate for the pronounced peak seen at n=60 in the reported mass spectra 
has been found. Apparently, the bond strain in Buckminsterfullerene is still more favorable than a 
situation with 20 unpaired electrons. Multiple-sheet models have not been considered in this work, 
since the number of dangling bonds in these would be even larger. 


Conclusions. 


We have demonstrated that new and unique information about large molecular systems can be 
achieved by using computational methods on a large scale. In many cases, that information would 
have been very difficult to obtain with other means. 

Even though our calculations on PAHs show a fairly slow convergence of properties, they allow 
for an extrapolation on the cohesive energy per carbon atom in bulk graphite. The data obtained 
from the Mulliken population analysis further show that the inner parts of the larger molecules 
studied have a domain with little total charge or spin density, which may mimic a graphite surface 
rather well. Since the problem of dealing with a large number of open shells can be circumvented 
by these calculations, the approach offers an attractive way of modeling a graphite surface as well 
as a way to study graphite intercalation compounds. 

The fact that pairing up dangling bond electrons may lead to more stable graphite fragments is 
somewhat surprising. This suggests that the edge of a graphite fragment may look different from 
what is commonly assumed. This observation may have far-reaching consequences for the under- 
standing of the chemistry and reactivity of large carbon clusters. 

Energetically, the lively debated soccerball-form of C 60 compares very favorably with a 
graphite-like sheet of carbon atoms. This, of course, does not prove that Buckminsterfullerene 
really is (fully or partly) responsible for the Cgo peak observed. Due to the large number of de- 
grees of freedom, a study of all plausible isomers would be a prohibitively cumbersome task. Even 
if such a complete coverage were possible the existence of different molecular species at these con- 
ditions is not likely to be governed by thermodynamic stabilities only. 
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Table 1 : Total energies and C-C bond distances computed for Cft n2He n and energies estimated 
using expression f2l with the parameters E c =-37.8363. E H =-Q.5685 Hartree. T he C-H 
distances were kept fix ed at 1.1 oA. 


n 

system 

' E tot,calc 

E tot,pred 

A 

Rcc (A) 

1 

^6^6 

230.4320 

230.4288 

0.0032 

1.384 

2 

C 24 H 12 

914.8938 

914.8932 

0.0006 

1.401 

3 

C 54 H 18 

2053.3942 

2053.3932 

0.0010 

1 .404 

4 

C 94 H 2 4 

3645.9003 

3645.9288 

-0.0285 

1 .406 

5 

C 150 H 30 

5692.4981 

5692.5000 

-0.0019 

1.406 a 


a assumed without geometry optimization. 
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Table 2: Total computed energies and distances for the diamond-like clusters, 
The C-H distances were kept fixed at 1.09A. 


n 

system 

-Etotfa.u.) 

Rcc (^) 

2 

C 10 H 16 

387.54332 

1.542 

3 

C 35 H 36 

1344.90665 

1.540 

4 

C 84 H 64 

3214.91903 

1.539 


Table 3: Computed properties for Si ycftm/nsterfa/terene, 


Bond lengths: 

R1 1.453 A 
R2 1.369 A 

Ionization potential: 

ASCF 7.92 eV 

HOMO energy 8.24 eV 

Electron affinity: 

ASCF 0.80 eV 

LUMO energy 0.60 eV 


Cohesive energy/carbon atom: 


114 kcal/mol 



